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Abstract 

The thermodynamic properties of heated nuclear matter are explored using 

an exactly solvable canonical ensemble model. This model reduces to the 

results of an ideal Fermi gas at low temperatures. At higher temperatures, the 

fragmentation of the nuclear matter into clusters of nucleons leads to features 

that resemble a Bose gas. Some parallels of this model with the phenomena of 

Bose condensation and with percolation phenomena are discussed. A simple 

expression for the hadronic equation of state is obtained from the model. 
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Properties of heated nuclear matter are of current interest, with heavy-ion collisions used 
to obtain information about nuclear matter away from the typical nuclear matter density 
and temperature. Theoretical concerns have centered around the nuclear equation of state, 
the specific heat, the behavior of the entropy per particle, and other thermodynamic is- 
sues ||T]-§|. The notion that nuclear matter undergoes a phase transition has gained special 
attention [|J. Most treatments of the equation of state have treated heated nuclear matter 
as a homogeneous system. Statistical p|-p3|] and percolation [p4j-|26|l models have been used 



by many groups to study nuclear multifragmentation and the role of clusters on nuclear 
properties. This paper explores various thermodynamic properties using an exactly solvable 
canonical ensemble model that allows for clusterization and leads to simple analytic expres- 
sions for thermodynamic quantities such as the equation of state. This canonical ensemble 
approach has properties that are similar in many ways to those obtained from models used 
in other areas, such as Bose condensation, Feynman's approach to the A transition in liquid 
helium and in polymer physics f2l| . For example, the cycle class decomposition of the 
symmetric group which appears in the symmetry of the bosonic system in Feynman's ap- 
proach is isomorphic to the cluster structure of a heated nuclear system. The cost function 
of moving a helium atom from one location to another is the same as internally exciting a 
cluster. 

The canonical ensemble partition function given below is obtained from the following 
simple picture. At T = 0, a system of A nucleons is in its ground state, which is treated as a 
degenerate Fermi gas. As T increases, the Fermi gas is excited into low lying excited states 
whose level density grows as p{E) ~ exp(2V aE). Particles and clusters can also be emitted 
with increasing T. The formation of clusters is an expression of the attractive nature of the 
nucleon-nucleon interaction. The repulsive nature of this interaction will be treated here in 
terms of a density dependent Skyrme approach At high density an excluded volume 

correction is needed to avoid some unphysical consequences of the model. At very large T 
clusterization is unlikely, and the nucleons behave essentially as a dilute Maxwell-Boltzmann 
gas. 



Considered as a statistical system, this model exhibits a phase transition in the fixed 
density but infinite A limit, at which the specific heat per nucleon is maximal. An infinite 
cluster suddenly appears at this point in a manner similar to the infinite cluster of percolation 
theory. Real nuclear systems are far from infinite collections of nucleons, but the evidence 
of such phase transitions can be obtained from a finite scaling analysis in the region of the 
phase transition [24], 4]]. 



In earlier papers [|20|-P3|1 the question of the thermodynamics of a fragmenting system 
was raised. There, it was assumed that the thermodynamic variables were contained in a 
single parameter x such that the weight of a particular cluster partition was proportional 
to x m , where m is the total number of clusters in that partition. Specifically, the partition 
function of the system was a polynomial in the variable x of degree A, where A is the number 
of particles in the system. This x contains physical quantities as discussed below in Eq. (^j). 
The computation of various thermodynamic quantities could then be reduced to functions 
of x and various moments of the number of clusters. Here the slightly more general case 
where there are two thermodynamic parameters x, y with the probability of a particular 
fragmentation outcome proportional to x m y A ~ m is analyzed. This choice reduces to the 
early models by noting x m y A ~ m = y A (x/y) m where x/y can be identified with the tuning 
parameter x used earlier and y A is an overall factor which has no effect on cluster yields. 
However, this more general case is necessary for an evaluation of some thermodynamic 
properties. 

We begin by assuming that each fragmentation outcome happens with a probability 
proportional to 

where x and y are functions of the thermodynamic variables V, T, (3k is the cluster size 
dependence of the weight discussed below, and n = (ni, n 2 , . . .) is the fragmentation vector, 
with rik the number of fragments with k nucleons such that J2k = A. The free energy 
for such a canonical system is given by F^(V,T) = — ksT In Z a with 



Z A = £ W(n) = £ Z^ ] 0)x m y A - m , (2) 

nen^ m=0 

where II^ is the set of all partitions of A and Z^™ is a function only of the vector (3, not of 
the thermodynamic variables. 

It is fairly straightforward to derive relations between logarithmic derivatives of the 
partition functions and moments of the multiplicity m = J2k n k, e.g. 

x-^\nZ A = (m), (x^j \nZ A = (m) 2 , ( 3 ) 

where {m)k = ((m — (m)) k ) is the fcth central moment. In general, (x-^) k In Z A {x) is the 
fc'th cumulant moment of m. The relations for derivatives with respect to y are similar. 

Given these formulae, the determination of thermodynamic quantities can be reduced 
to a computation of the cumulant moments of the fragmentation multiplicity m, since only 
logarithmic derivatives of the partition function (such as those given above) are used in 
the determination of the energy, pressure, specific heat, etc. Applying the usual relations 
between the thermodynamic functions and the partition function, the result for the specific 
heat and pressure is 

CV_ x (^^_ (TdxV T 2 d 2 x 
~k^~^[~df \~x~dT) + ~df 2 



+ (A - (m)) 



2T dy (TdyY T 2 d 2 y 



y dT \ydTJ y dT 2 

+ ((m 2 ) - (m) 2 ) 



'Tdx_ _ Tdy_\ 2 
~x&T ~ ~y~&f) 



p y , ,V dx , A . „V dy 

where the last equation is the equation of state. The calculation of the cumulant moments 
of m can be done by applying Eq. (|3|) and using a recursive evaluation of the coefficients of 
the partition function Z A m \ namely |23| , 

4 m) = -j-E , (5) 
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where Z A = 1//3a- The whole recursive procedure is easy to implement as a computer 
program. 

The choice of x, y and (3k is determined by the physics of the situation. The term x m 
comes from phase space factors and translational partition function considerations, namely 
x = V/\ d , with A the thermal wavelength of a nucleon and d the dimensionality of the 
system (with d = 3 appropriate for nuclear fragmentation). This x, with (3k = k l+d l 2 
in Eq. ([!]) is essentially the weight given to Bose condensation problems in d dimensions. 
In this case, n& is the number of cycles of length A; in a cycle class decomposition of the 
permutation associated with a particular Bose gas state. For a nuclear system, the term 
y is due to cluster binding and internal excitations. In a simplified view of binding energy 
considerations, each fragment of size k has a binding energy of ay(k — 1), so a total binding 
energy of Eb ~ ay (A — (m)) suggests a Boltzmann weight of y — exp(av/fcsT). It is 
appropriate to also include in this factor a Fermi gas level density term arising from internal 
excitations, leading to the result, 



where ay is the binding energy per nucleon, Eq is the level spacing parameter for excited 
states, and T is a temperature cutoff factor for internal excitations. The expression x/y is 
the tuning parameter in previous papers p0|-p3]]. 

In general the parameters ay and 6q are density dependent. In the Skyrme approach, 
ay(p) is given by 



and for a Fermi gas, Eo(p) = (A/tt 2 )ef(p) with ef(p) ~ p 2 ^ 3 the Fermi energy, a^ = 3ef(po)/5 
and ao, 03 are Skyrme parameters, which can be determined by requiring the binding en- 
ergy/nucleon to be a maximum at p = po with the value of 8.0 MeV. It should be noted that 
the factor exp(/cgT / e ) can be rewritten as exp(a(md 2 /2/2, 2 )fc^T) . Using the above relation 




(6) 
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between Fermi energy and the density of states, a is a numerical constant close to 1 and 
d 3 = V/A. This factor with a = 1 is the cost function of moving a Helium atom from one 



location to another in Feynman's approach to the A transition in liquid Helium |27| . 

The only remaining parameter to set in this model is (3k- If the weight W(n) from Eq. [I] 
was for a Bose gas, then |3B| 

Xk = ¥ + \' (8) 

with r = 1 + d/2. The 1/k term arises from the zero momentum states of the Bose con- 
densate. This additional term is irrelevant at high temperatures, but of some importance 
below the critical point. Since the low temperature behavior we want is essentially Fermi 
gas like, we can ignore the zero momentum term, which suggests we use (3^ = k T . Here r is 



the critical parameter r discussed in Fisher [31|, Finn, et. al. [32] and Gilkes, et. al. B. We 



choose r = 5/2 to match the exponent of a Bose gas in three dimensions, but data from 0] 
suggests a somewhat lower r w 2.2. Although in a Bose gas r is fixed by the dimension of 
the system, this nuclear fragmentation model is not so constrained, and a different r can be 
used. 

Applying the above choice of parameters to Eq. ([|) gives 

C v , . d , A , * \ 2k B T ( T ^ 3 



/ \ d , (A I \ \ B ( 1 \ 

k B {m) 2 + {A - {m)) —[TTTj 

(d a v k B T ( T \ 
PV = (m)k B T + (A- (m)) 



2 , s 5tt 2 / k B T \ 2 T 



x -e F (p) 1 + 



5 yrJ \ 12 \e F (p)J T + T 0/ 

_ ao Q +as(1 + ff) Q ,+ °). (9) 

Note that the equation of state is particularly simple. The first term is the ideal gas law, 
while the second term contains the low temperature degeneracy pressure and interaction 
terms from the Skyrme potential. 
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The specific heat and equation of state given by the above expression are plotted in 
Fig. |1|. For low T, Cy oc T and is that of a heated nucleus of nucleons treated as a nearly 
degenerate Fermi gas. At higher T, nucleons and clusters are emitted from this nucleus, 
increasing the number of degrees of freedom and the specific heat. At a relatively low 
temperature this effect causes a transition from a degenerate Fermi gas to a Bose gas-like 
state. This change in Cy shows up in the figure as a shoulder. At very high T the nucleus 
of A nucleons dissolves into A independent nucleons and Cy = \dksA. The "cusp" like 
behavior of Cy can be understood as a critical point in the infinite A limit as discussed in 
the next paragraph. A rounded peak is seen instead of the cusp because of the finite size of 
the system. Figure |I](b) illustrates the behavior of P with V/A in the transition region. In 
this region, the equation of state behaves in many ways like a van de Waals gas. 

The thermodynamic limit A, V — > oo with p = A/V finite is of interest since it specifies 
the critical point behavior and the nature of any phase transitions. In this case, we can work 
in the grand canonical limit and consider the question of behavior of the largest cluster. A 
specific concern is whether this theory has the same characteristic as the infinite cluster in 
percolation theory i.e. for p > p c the infinite cluster exists while for p < p c it does not 



exist. In the grand canonical limit (n^) is given by 



= y— = yW (10) 



which implies the mass constraint 

A 



Xyjfc 1 -' (11) 



x /y k 

Here Aj '(x/y) is finite in the thermodynamic limit since x oc V and y = y(p,T). For 
z < 1, the sum is always less than or equal to the case z — 1, for which the sum gives 
£(t — 1), i.e. the Riemann zeta function. When z > 1 the sum diverges and therefore the 
constraint cannot hold. Now Aj '(x/y) is a function of A/V,T which are fixed, so z is the 
only parameter we can adjust. If (x/y) /A < 1/((t — 1) then we can find a z < 1 such that 
the mass constraint is satisfied. However, if (x/y) /A > 1/C( r — 1) then the constraint can't 



be met. The sum must be truncated, which implies that the expectation of certain large 
clusters must be identically zero. In other words, when (x/y)/A > 1/C( r — 1) 5 there can be 
no infinite cluster, and (x/y) c /A= l/£(r — 1) defines a critical point for this system, which 
is identical to the critical point of a Bose gas if y — 1. Since ((r — 1) < oo only if r > 2, 
this also implies that the infinite cluster can not exist if r < 2. 

Universality implies the model near the critical point is specified uniquely by two critical 
exponents. We already mentioned the critical exponent r, but another exponent is needed 
to determine the critical behavior. By the above discussion, the fraction of mass not in the 
infinite cluster, m x = lim^oo 12£=i k{nk)/A satisfies m x < 1 for T < Tc and m x = 1 for 
T > Tc- Thus (rioo) = 1 — m x is zero for T > T c . Near T = T c , (tt-oo) ~ (T — Tc)/Tc- 
The analog of m x in the Bose system is the fraction of Bose particles in excited states, while 
(rioo) is the fraction in the condensed ground state. For the Bose system the order parameter 
is taken as the square root of the number of Bosons in the ground state, and in analogy 
\J (rioo) oc (T — T^) 1 / 2 , which gives the second critical exponent (3 = 1/2. 

In summary, this paper introduced a simplified model for nuclear systems which attempts 
to be valid across a wide range of temperatures and densities. Building on a model devel- 
oped earlier |2Ti|-|2"5[ with simple expressions for cluster yields, fluctuations and correlations 
useful in the analysis of inclusive and exclusive data, the thermodynamic properties of the 
model have been made explicit and the corresponding thermodynamic functions such as 
the pressure and specific heat were shown to be simple analytic functions of the density, 
temperature and cumulant moments of the multiplicity. At typical nuclear matter densities 
and very low temperatures the model is equivalent to a nearly degenerate Fermi gas. At 
very high temperatures and low densities it reduces to an ideal Maxwell-Boltzmann gas of 
nucleons. In intermediate regions, the nucleus will break up into clusters of various sizes and 
the system has some features similar to that of a Bose gas. More importantly, this model 
exhibits a phase transition in the infinite A but finite density limit with critical parame- 
ters similar to other models of nuclear fragmentation. As in percolation theory, the infinite 
cluster appears suddenly at the critical point. As in a Bose gas or Feynman's model of the 
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A transition in liquid helium, the specific heat has a maximum at the critical point with a 
discontinuous derivative. The "cusp" in this model however is smoothed out by the finite 
size of the system. 

This work supported in part by the National Science Foundation Grant No. NSFPHY 
92-12016. 



9 



REFERENCES 

[1] G. Bertsch and S. DasGupta, Phys. Rep. 160, 189 (1988). 
[2] L. Csernai and J. Kapusta, Phys. Rep. 131, 223 (1986). 
[3] S. DasGupta and G. Westfall, Phys. Today 46(5), 34 (1993). 
[4] M. L. Gilkes, et. al., Phys. Rev. Lett. 73, 1590 (1994). 

[5] H. R. Jaqaman, A. R. DeAngelis, A. Ecker and D. H. E. Gross, Nucl. Phys. A541, 492 
(1992). 

[6] H. R. Jaqaman and D. H. E. Gross, Nucl. Phys. A524, 321 (1991). 

[7] A. R. DeAngelis and A. Z. Mekjian, Phys. Rev. C 40, 105 (1989). 

[8] A. R. DeAngelis, D. H. E. Gross and R. Heck, Nucl. Phys A537, 606 (1992). 

[9] D. H. E. Gross, A. Ecker and A. R. DeAngelis, Ann. d. Phys. 1, 340 (1992). 

[10] D. H. E. Gross, L. Sapathy, M. Ta-Chung and M. Sapathy, Z. Phys. A309, 41 (1982), 
D. H. E. Gross and X. Zhang, Phys. Lett. 161B, 47 (1985). 

[11] S. Sobotka and L. Moretto, Phys. Rev. C31, 668 (1985). 

[12] J. Aichelin and J. Hufner, Phys. Lett. 136B, (1984). 

[13] S. Koonin and J. Randrup, Nucl. Phys. A474, 173 (1987). 

[14] D. H. Boal and J. N. Glosi, Phys. Rev. C 38, 1870 (1988). 

[15] A. J. Cole, Phys. Rev. C 40, 2024 (1989). 

[16] C. Fai and J. Randrup, Nucl. Phys. A404, 551 (1983). 

[17] W. A. Friedman, Phys. Rev. C 42, 667 (1990). 

[18] J. Bondorf, R. Donangelo, I. N. Mishustin, and H. Schule, Nucl. Phys. A444, 460 (1985). 



10 



[19] K. Sneppen, Nucl. Phys. A470, 213 (1987). 

[20] A. Z. Mekjian, Phys. Rev. Lett. 64, 2125 (1990); Phys. Rev. C 41, 2103 (1990). 
[21] A. Z. Mekjian and S. J. Lee, Phys. Rev. A 44, 6294 (1991). 

[22] S. J. Lee and A. Z. Mekjian, Phys. Lett. A149, 7 (1990); Phys. Rev. C 45, 365 (1992); 
45, 1284 (1992). 

[23] K. C. Chase and A. Z. Mekjian, Phys. Rev. C 49, 2164 (1994). 
[24] X. Campi, J. Phys. A19, L917 (1986); Phys. Lett. B, 351 (1988). 
[25] J. Desbois, Nucl. Phys. A466, 724 (1987) 

[26] W. Bauer, D. R. Dean, U. Mosel and U. Post, Nucl. Phys. A452, 699 (1986). 

[27] R. P. Feynman, Statistical Mechanics: A Set of Lectures, Frontiers in Physics (Ben- 
jamin/Cummings, Reading, MA, 1972). 

[28] F. P. Kelly, Reversibility and Stochastic networks (Wiley, New York, 1989). 

[29] M. W. Cortin, H. Toki and D. K. Scott, Phys. Lett. 123B, 289 (1983). 

[30] H. Jaqaman, A. Z. Mekjian and L. Zamick, Phys. Rev. C 27, 2782 (1983). 

[31] M. E. Fisher, Physics (N.Y.) 3, 255 (1967). 

[32] J. E. Finn et. al, Phys. Rev. Lett. 49, 1321 (1982). 

[33] K. Huang, Statistical Mechanics (2nd ed.), (Wiley, New York, 1987). 

[34] D. Stauffer and A. Aharony, Introduction to percolation theory (Taylor and Francis, 
London, 1992). 



11 



FIGURES 

FIG. 1. The specific heat (a) and equation of state(b). The specific heat is Bose-like except at 
low temperatures, where it transitions to Fermi gas like behavior. In this region, the equation of 
state shows the characteristic van de Waals behavior. 
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